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A major difficulty in comparing quantum and classical behavior resides in the struc- 
tural differences between the corresponding mathematical languages. The Heisenberg 
equations of motion are operator equations only formally identical to the classical equa- 
tions of motion. By taking the expectation of these equations the well known Ehrenfest 
theorem provides identities which, however, are not a closed system of equations which 
allows to evaluate the time evolution of the system. The formalism of the effective action 
seems to offer a possibility of comparing quantum and classical evolutions in a system- 
atic and logically consistent way by naturally providing approximation schemes for the 
expectations of the coordinates which at the zeroth order coincide with the classical 
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evolution [1]. 

The effective action formalism leads to equations of motion which differ from the 
classical equations by the addition of terms nonlocal in the time variable. This means 
that for these equations an initial value problem is not meaningful and they have to 
be interpreted in an appropriate way. Here we analyze situations in which the nonlocal 
terms can be reasonably approximated by local ones so that the quantum corrections do 
not modify the locality of classical equations. In the simplest approximation, the effective 
Lagrangian differs from the corresponding classical one by a renormalization of both the 
potential and the kinetic energy terms. We shall not discuss the causal formalism used, 
for example, in Refs. [2-4], as in the approximation considered this would lead to the 
same local equations. 

The present contribution describes the beginning of a systematic study of semiclassical 
evolutions using the effective action formalism. In the first part, after introducing the 
formalism of the effective action and its expansion in powers of h (loop-expansion) in 
the context of quantum mechanics, we concentrate on the structure of the first order 
corrections in ft. These corrections are evaluated to the second order in the derivative 
expansion [5], by two different methods. The first is based on a Euclidean approach [6], 
the second one on an adiabatic approximation in evaluating functional determinants. 

In the second part of the article we put the formalism at work, choosing as our case 
study a two-dimensional (2-D) anharmonic oscillator of the kind considered in molecular 
physics. The results of the simulations show that by increasing ft the effective dynamics 
tends to regularize the classical motion and becomes qualitatively very similar to the 
quantum evolution provided the energy is sufficiently small. 

The evaluation of the effective dynamics in more general cases will be presented in a 
forthcoming paper. 

PART I 

1. Effective action in quantum mechanics 

In this Section we define the effective action [7] . For simplicity, consider a one degree 
of freedom Hamiltonian, 

(1.1) H{p,q)^Ho{P,q) + U{q), 
where 

(1.2) H^(^p^q)^I_ + lmuj^^ 

2m. 2 

and the confining potential U{q) is an even polynomial of q. We choose the constant of 
U{q) so that the lowest eigenvalue of H is 0. The generating functional of the Green 
functions is 

(1.3) z[j] = (o|r(e* ^ dtj{t)m)\Q) 
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where |0) is the ground state of iJ, q{t) — e^^''^ qe^^^'^ , J{t) is a source vanishing for 
\t\ -^ oo and T is the time-ordering operator. In Eq. (1.3), as well as in the following, 
the integrations with boundaries not explicitly indicated are to be understood between 
— oo and +00. The generating functional of the connected Green functions is defined as 
W[J] = ~ihlnZ[J] and the Legendre transform of W[J] gives the effective action. By 
indicating with q the variable conjugated to J, i.e., 

we define 

(1.5) T[q] = W[J]- Jdtq{t)J{t), 

where J has to be thought, inverting relation (1.4), as a functional of q. The functional 
r[q] represents the analog of the classical action, S[q] ~ J dt (^^■mq^{t) — V{q{t))'^, where 
V{q) = ^muP'q^ + U{q), and can be written in the form 

(1.6) V[q] = S[q]+fn[q], 

with T{)[q] = 0. The Legendre transform can be calculated using the methods of Ref. [8]. 
From the functional derivative of the classical action with respect to the position q{t) 
one obtains the Euler-Lagrange equation of motion 

<-' H - -<"^ 

In the same way the functional derivative of the effective action T[q] with respect to the 
q{t) given by (1.4) yields 

This equation can be rewritten in the form 

(1.9) ^q^t) + d,V{q{t))-^-^ = J{t). 

As we shall see in the next Section, Th\q\ admits an expansion in powers of fi, whose 
coefficients have a simple diagrammatic interpretation (loop-expansion). In this way we 
can view the quantum integro-differential equation (1.9) as a perturbation of the classical 
equation of motion. 

In order to interpret the solutions of Eq. (1.8), we rewrite Z[J], defined in (1.3), in 
the equivalent form 

(1.10) Z[J] = (0|t//(+oo,-oo)|0) = (0|r(e-^/'**[^-«^(*)l)|0), 



4 Fabrizio Cametti, Giovanni Jona-Lasinio, Carlo Presilla and Fabio Toninelli 

where Ug{ti,, to) is the evolution operator from ta to % in the Schrodingcr representation 
and in presence of the external source J{t). Note that in quantum mechanics J{t) is an 
external force. The variable conjugated to J then is 



("[jjIAj]) 



where 



(1.12) (a[,,]| = (0|C/i(+cx3,0) 

(1.13) |/3[j]) = [/i(0,-cx3)|0). 

Since generally |/?[j]) differs from |q;[j]), q{t) is a nondiagonal matrix element of q between 
two states which evolve in presence of Jit). The solution of Eq. (1.8) therefore can be 
complex valued. In the harmonic case U{q) constant, \oL\n) and |/3[ji) are coherent states 
and coincide, up to a phase, if J{uj) ~ 0, where J is the Fourier transform of the external 
force. In the anharmonic case, more complicated conditions have to be imposed on J 
so that the two states coincide. If these conditions are satisfied, q{t) is the expectation 
value of the position operator. 

2. Loop expansion of the effective action 

The effective action cannot be evaluated exactly for anharmonic systems, i.e., U{q) ^ 
constant. A widely used approximation scheme is the loop expansion (see for example 
[9,10]), or semiclassical approximation, consisting in an expansion of T[q] in powers of h. 
At the lowest order the effective action coincides with the classical action, whereas the 
one-loop term is expressed by means of a functional determinant. 

In order to obtain the loop expansion we express Z[J] as a path integral. Equation 
(1.10) can be rewritten, using the Gell-Mann and Low theorem [11], as 

(2.1) Z[J] = >"' ^ ^' °' 



(Oo[r(e-i/'^*^(9o(*)))|0o) 



where (7o(i) = e'i^°*qe i-^"* and [Oq) is the ground state of Hq- Equation (2.1) is 
equivalent to 



(2.2) Z[J] 



e-i/Ml7%)zo[J'] 


J'=J 


e ^ -^ "^'^l "'(*)) Zo[ J'] 


J'=0 



where Zo[J] in terms of Feynman path integrals [12] reads 

(2.3) Zo[J]= lim /dxdy /%]^e^-f-^''*[^«'(*)-^"'«'(*)+'^(*)«(*)Vo(:c)^o(y). 



T-+00 
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Here (/3o(2;) = (2;|0o) and d[q\y is the functional measure on paths with endpoints q{—T) = 
X, q{T) = y. The purely oscillating integrand in Eq. (2.3) can be regularized by changing 
Lu into ujg = uj{l — ie) with e ^ 0+ [13]. Comparing Eqs. (2.2) and (2.3) we obtain 

(2.4) Z[J] = lim lim 

e^O+T^oo 

JdxdyJ d[q]v e^ ^-T dt[^e(t)-^u.U\t)-u(,(t))+j(tUt)] ^p(^)y,„(^) 

/dxciy/d[g]^e^'''-^''*[^«'(*)-^"?«'W-^(«(*»]^o(:c)^o(2/) 

Now we apply the stationary phase approximation to (2.4), expanding the exponent 
at the numerator around the solution qoit) of 

(2.5) mqa{t) = -muj%{t) - dgU{qo{t)) + J{t) 
which vanishes for \t\ -^ oo. We find 

(2.6) Z\J] ~ e"s^"^[®l+''^''*'^(*^*'^*^) 

X lim lim 

e^O+T^oo r£;[g]0g^/5T''*[9^(*)-'^?9^(t)] 

Note that the integrations over x and y disappear since ipo{-) is proportional to S{.) in 
the limit /i — > 0. The Gaussian integrals in (2.6) can be performed yielding 

(2.7) Z[J] ~ e^(-^[*l+'''''*-'(*^*(*^) 

e^0+ T^oo I det {-d^ - ujf) I 

where the differential operators act on functions y(t) with Dirichlet boundary conditions 
y{-T) = y{T) = 0. From Eq. (2.7) we obtain 

(2.8) W[J]^Wo[J] + hWi[J]+0{h^) 

^S[qo\+ I dtJ{t)qo{t) 



5^0+ T^oo 2 \ det {-dt - tj|) / 



Setting q = qo + hq' and remembering that -r-jrl 
one-loop order is 

(2.9) T[q]^To[q] + flT,[q]+0(fl^) 

= Wo [J] + hWi [J] - I dtq{t)J{t) + O {ti 



— —J{t), the effective action to 

90 
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S[q - hq'] + hWi[J] - h I dtq'{t)J{t) + O (h^) 



2 e^o+ T^oo \ det (-9j - w|) / 



-o(r). 

■ichlct(±T) 



3. Derivative expansion of the effective action 

The classical action S[q] is the time integral of a density (the Lagrangian) which is 
an ordinary function of q{t) and q{t). As a consequence, the classical equation of motion 
(1.7) is a differential equation. On the other hand, the effective action T[q] is nonlocal 
in time and, therefore, the variational equation (1.8) is also nonlocal. If q{t) varies 
slowly, however, it is possible to expand T[q] around a constant value of q (derivative 
expansion [6, 14]). In this expansion one finds that also T[q] can be written as the time 
integral of a density, which is a series of terms involving time derivatives of q{t) of 
increasing order: 

(3.1) r[q] =Jdt (^-VMt)) + ^^^e{t) + A{q{t))qHt) + B{q{t))if(t) + . 

As we shall see, the derivative expansion (3.1) does not generally converge and has only 
an asymptotic validity for q{t) -^ constant. The absence in (3.1) of odd powers of q{t) is 
a consequence of the time reversal symmetry of the Haniiltonian. 

Except for Ve and Z, all the terms in the derivative expansion (3.1) are at least of 
order h: 

(3.2) VM = Irnq'^LU^ + U{q) + hV,i (q) + O (ft^) 

(3.3) Z{q) = m + hZi{q) + (h^) 

(3.4) A(q)=nAi{q) + 0(h^) 

(3.5) B{q) = hB^{q) + 0{h^). 

The effective potential 14(9), well known in quantum field theory in the study of sponta- 
neous symmetry breaking [5], is everywhere convex [15]. It may happen that the effective 
potential evaluated at a finite h order loses somewhere its convexity if the classical po- 
tential is not everywhere convex [16]. In this paper we restrict ourselves to a phase-space 
region where the evaluated effective potential is convex. 

If the derivative expansion (3.1) is truncated at a finite order 2iV, the corresponding 
variational equation is a differential equation of order 27V. We thus have a Cauchy 
problem with 2N initial conditions. It is clear that these conditions do not determine 
completely the initial wave function of the system. They are constraints which must be 
imposed in the choice of the initial wave function for a comparison between true and 
effective quantum evolutions. We confine ourselves to the second order in the derivative 
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expansion (DE2), that is 



(3.6) 



TM 



dt\-vMt)) + ^^e{t) 



This is the simplest approximation to the effective action which preserves the structure 
of the classical equations of motion. 

In the following we work out and compare two methods to obtain the derivative 
expansion of the effective action. The first is an adaptation to quantum mechanics of a 
method [6] used in quantum field theory and based on the Euclidean functional formalism. 
In the second method, we relate the derivative expansion to the adiabatic approximation 
of a differential equation with slowly varying coefhcients. In this way we are able to give 
an estimate of the validity of the derivative expansion. 

3'1. Derivative expansion: Euclidean approach. - The derivative expansion of the 
effective action can be obtained starting from the Euclidean generating functional 



(3.7) 



Ze[J] = 



/d[g]|e-*(^^M-.H*''(*)«W)(po(x)¥'o(2/)da;dj/ 



where the Euclidean action ^^[g] is defined by 



(3.8) 



SE[q] 



dt ( ^mq\t) + ^mcuWt) + U{q{t)) 



Setting VK_b[J] = ft In Z^; [J] and q{t) = T^Kt)^ ^'^ introduce the Euchdean effective 
action 

(3.9) TE[q] = We[J] - I dtJ{t)q{t). 

In analogy with the results of Section 2, to one- loop order we have 



(3.10) 



^E[q] = -SE[q] - 77 In 



det 


r <5'Se[9] 1 


Sq(t)8q(s) 


det 


■52 5e[0] 


[Sq(t)Sq(8)\ 



-Oih'), 



where the differential operator ^ (t)s(s) ^^^ ^^ rewritten as 



(3.11) 



S'SeM 



= [-md^ + mu;^ + d^U{q{t))] 5{t - s). 



5q{t)5q{s) 
The second order in the derivative expansion of the Euclidean effective action is 



(3.12) 



^E[q]^- j 



dt 



VMt)) + \z{q{t))q\t) 
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where Ve{q) and Z{q) are the same functions that appear in (3.1). 

The effective potential can be found by combining (3.10) and (3.12) for q{t) constant 



(3.13) 



I dtVei{q)^^lndet 



S^Ssiq] 



Sq{t)Sq{s) 



In det 

2 



Sq(t)6q{s) 



Employing the functional analogue of the identity In det A = tr In A, valid for any Her- 
niitian matrix A, we get 



(3.14) 



In dot 



S'SeM 



We use the Dirac notation to write 



tr In 



S^SeM 



6qit)5q{.s) 



(3.15) 



1 



d^ + Lu' + -d'U{q) 5{t - s) = {t\ \P'+u' + -d'U{q) \s), 



where the P operator is defined in the {\t)} basis by {t\P\s) — —i-^5{t — s). Equation 
(3.13) becomes 

(3.16) j dtVei {q)^^fdt( {t\ In (p^ +u;^ + l^d^gU{q)] \t) - {t\ In (p^ + cj^^ \t)\ . 
With the help of the identity J dp\p){p\ = 1, where P\p) ~ p\p), we can write 



(3.17) (i|ln [P'+uj' + -d',U{q)) \t) = ± fdp\n (p'+oj^ + -d'^U(q) 



The integral in the above expression can be evaluated exactly and from (3.16) we finally 
get 



(3.18) 



V^M = \[\h^^ + yiU{q)-u: 



The determination of Zi{q) is more involved. From Eq. (3.12) we see that Z(c[) is 
the coefficient of the term containing (7(i)^ in the effective action. We can thus write 



(3.19) / dtZi{q{t))q\t) = Indct 



5^SE[q] 



5q{t)5q{s) 



- In det 



5^SE[qc 



5q{t)5q{s) 



qa^q(t) 



with the assumption that we consider in the r.h.s. only those terms with at most two 
time derivatives of q{t). The first term in the r.h.s. of (3.19) is essentially the one-loop 
Euclidean effective action, while the second one comes from the effective potential. The 
second functional determinant in (3.19) has to be evaluated with a constant qc which, at 
the end, must be replaced with q{t). The terms due to the normalization of Ze[J], being 
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common to both the effective action and the effective potential, cancel each other. Again 
we change the logarithm of the determinant into the trace of the logarithm and write the 
differential operators in Dirac notation. It is useful to introduce two operators, P and T, 
satisfying the commutation relation [T, P] = i, and with elements (i|P|s) = —i-^5{t — s), 
{p\T\q) = i-§-;S{p — q), where \t),\s) and \p), \q) are eigenstates of T and P, respectively. 
In addition, we write the difference of two logarithms of positive defined operators in the 
parametric form 

(3.20) \nA-\nB= r^ (g-^^ - e'^' 
We then arrive at the following expression 

(3.21) fdtZMt))q^t) 
dtlt\ n ^e"(™-^'+""'+^''^(«(*)))" - g-(mP^+mc^^+a,^C/(<j(f)))s| A 

Jo s ' 
Since we keep only terms at most quadratic in qit), we can expand d'^U{q{T)) as follows: 

(3.22) dlU{q{f)) - dlU{q{t)) + Qait) + \QH{t), 

where Q = f-t, a{t) = dtd'^U{q{t)) and b{t) == df d'^U {q{t)) . All the terms proportional 
to Q", with n > 3, are neglected since they do not contribute to the determination of 
Zi{q). The expressions for (ile-^-^'^li) and for (i|e-['"-P'+<5'^(*)+5Q'Kt)]«|i) ^re known 
[17] and can be inserted in (3.21). Finally we can expand the integrand in Eq. (3.21) 
maintaining only the terms linear in b(t) and at most quadratic in a{t). Performing the 
integration over the variable s, we obtain 

(3.23) ZM ' <^?^<'»^ 



32m2 



{u;^ + ±d^^U{q)y 



3'2. Derivative expansion as a WKB-like approximation. - The functional determinant 
in the one-loop term of the effective action (2.9) can be expressed by means of the 
Gelfand-Yaglom formula [18-20] as 

(3.24) r[q] = S[q] + ^ lim lim In ( '^f'^'^^, ) + O (h^) 
where Fg (t) is the solution of 

F,(i) + K + ^a2c/(9(t)))F,(t)^o 

(3.25) <( F,i-T)=0 

M-T) = 1. 
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Note that the time variable appearing in the above equations is the real time. At first 
sight it might seem that, on account of the factor i in (3.24), the one-loop contribution to 
T[q] is imaginary if q{t) is real. Actually, the effect of the rcgularization a; ^ oj^ is such 
that ri[g] has generally both a real and an imaginary part. As we shall see, the latter 
disappears if q{t) varies slowly with time. Without the regularization the expression for 
ri[g] would be ill-defined, both numerator and denominator oscillating with T. 

We obtain the derivative expansion of the effective action at order h starting from Eq. 
(3.24). For the moment we neglect the frequency regularization which we will reintroduce 
later. In order to deal with convergent integrals we suppose that q(t) — for \t\ > s. At 
the end of the calculation, i.e., after the limits T ^ oo and e ^ 0+ have been taken, we 
will let s — ^ oo. 

If we set q{t) = Q{pt), q{t) varies slowly if p is small. The expansion of r[(7] around q{t) 
constant is therefore related to the asymptotic expansion of F{t) for p ^ 0. Introducing 
the variable t = pt and setting <I>(t) = F(t/ p) and A:^(r) = w^ + (1/m) d'^U(Q(T)), Eq. 
(3.25) becomes 



(3.26) 



^^r) + j,e{T)^{T) 


= 


<^{-pT) = 




1-r'^i-pT) - i. 





An approximate solution of (3.26) for p ^ can be found by means of the WKB method 
[21] with the parameter p playing the role of ft. The A'^-th order solution is 



(3.27) $,«(t) 



1 



v^w^ 



c+e 



J /" W2n(T')dT' ^ c_^--,.r W'2N(r')dr' 



where W2n(t) is obtained, neglecting all the terms of order higher than p^^, from the 
recursive relation 



(3.28) W^^ir) 



k\T)+p'^W2iN-l){T) 



dT^ \^/W2{N-1){t) 



with 



(3.29) 



Wo{t) = fc(T). 



Imposing the initial conditions and going back to the variable t we find that at the lowest 
order the solution of (3.25) is 



(3.30) 



Foit) = 



^tf'_^dtyu^^ + ^d^U{q{t')) _ ^tJ'_^dt'y/u^^ + J^dlU(qit')) 



2i^uj^uj^ + ^dp{q{t)) 
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When the regularization a; — > Wg is reintroduced and -Fo(0 i^ evaluated at the time T, 
the second exponential, proportional to ^^'^'^^At-s) ^ vanishes for large T and can be 
neglected, since the limit T ^ oo has to be performed before the limit e ^ 0+. We 
obtain therefore 



(3.31) T[q] ^ S[q] ^ f^ dt i ^c.2 + ^dlU{q{t)) - c j . 

Recalling (3.2), the first quantum correction to the classical potential is 



1 



(3.32) v,M = ^\\j^^ + -dlU{q)-u: 

The next order of the WKB approximation gives 




(3.33) T[q] c, S[q] -■- I dt\ Juj^ + -a2C/(g(t)) - .. 



m 



2 



2 / /*3^— — -'^ (*) 



vi^hich implies 

(3.34) Zi{q) - 



{cu^ + ±dmiqit))y 



{d^M{q)f 



Equations (3.32) and (3.34) agree with the results found in Section 3T. 

If the classical potential V{q) is not everywhere convex, in the regions where uj"^ + 
^dgU{q) is negative the effective potential Vei{q) and Zi{q) become imaginary. More- 
over, Zi(q) has a divergence at the points where where uj'^ + ^d'iU{q) = and this 
corresponds to the fact that the WKB approximation loses its validity near the turning 
points fc^(T) ~ 0. 

It is clear that the iV-th order WKB approximation for $(t) corresponds to the 
derivative expansion of Ti[q] at order 2N. One can also check that no terms with an odd 
number of derivatives appear. The connection to the WKB approximation also shows, 
as previously stated, that the derivative expansion has only an asymptotic validity for 
q^O. 

From Eq. (3.24) it is clear that iiq{t) is real, r[q] up to one-loop order is not necessarily 
real. However, from Eq. (3.28) we see that, if the classical potential is everywhere convex, 
all the terms of the derivative expansion of the effective action are real if q(t) is real. The 
contradiction is only apparent. It can be seen that the imaginary part of T[q] is due to 
singularities in the Green functions which do not contribute to the derivative expansion. 
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In conclusion, in the case of vanishing external source J(t) = the DE2 approximation 
at order h of Eq. (1.8) reads 

(3.35) (m + hZ,{q{m m + ^dqZMm'it) = -dq {V{q{t)) + hV,^{q{t))) . 

We discuss the validity of this equation in the case V{q) — ^muo'^q^ + j^q'^- Equation 
(3.35) is approximate both because the DE2 approximation is adopted and because the 
terms of order higher than h are neglected. For a solution q{t) of amplitude A these two 
approximations are valid if 

(3-36) i'' 242 « 1 

and 

(3.37) -^ < 1, 

respectively. Under these conditions the solutions of Eqs. (3.35) and (1.8) remain close 
for a time t satisfying 

(3.38) ci 4^ ^M_ « 1. 

PART II 

4. 2-D anharmonic oscillator: classical 

Classical systems with more than one degree of freedom present a richer variety of 
phenomena and in particular they may exhibit chaotic behavior for J — 0. The formalism 
described in Part I can be generalized without difhculties to many degrees of freedom. 
Here we study the system whose Lagrangian is [22] 

(4.1) L{qi,q2,qi,q2) == -^miql + qj) - -muj'^{qf + ql) - gqlql- 

Apparently, the system has four free parameters: m, lu, g and the energy E. However, 
the rescaling t -^ t/uj, qi -^ qi^ muP' j g , qi -^ qi^moj^/g, for i ~ 1,2, yields 



(4.2) L 



? 4 

m LU 



1 / -2 I •2\ -'- / 2 I 2\ 2 2 

2(91 +12)- 2(91 +<i2)-qiQ2 



where, now, qi, qi and t are dimensionless. The energy of the system (4.1) is then 
E ~ {m'^uj*/g)s, where s is the dimensionless energy of the dimensionless Lagrangian 
-^ = ^{Qi +92) ~ 5(91 +92) ^ Qi^li- We conclude that e is the unique free parameter of 
the system under consideration. 
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The rescaled equations of motion 



(4.3) 
(4.4) 



91 

92 



-qiil + lql) 
-92(1 + 2g2) 



have been numerically integrated using a standard fourth order Rungc-Kutta method 
[23] . A qualitative description of the corresponding solutions has been achieved by con- 
structing the surfaces of section (Poincare sections) [24] and evaluating the largest Lya- 
punov exponent [25] . The degree of chaoticity of the system can be summarized by the 
fraction of regular orbits on the energy shell as a function of the dimensionless energy 
e. This fraction is close to unity for e < 0.75 and vanishes exponentially for e > 0.75. 
The border value e ~ 0.75 agrees with that obtained from the Toda criterion [26]. In 
our system, the sign of the curvature of the energy surface where the motion takes place 
is given by sign(det He(y)), where He(y) is the Hessian of F = {q\ + 9i)/2 + g^gl, 
i.e., He(y)ij = J^ Y . ■ This sign changes from positive to negative at e == 3/4. It is 



worth noting that the Toda criterion does not detect the first occurrence of chaos in 
the Poincare sections [27,28]. Nevertheless, when e > 3/4 we find that chaotic orbits 
are spread all over the sections. For e < 3/4, the irregular orbits are located in a small 
region of the Poincare sections, namely near a perturbed separatrix where chaos initially 
appears in consequence of the mechanism of the heteroclinic intersection [29] . 



5. 2-D anharmonic oscillator: quantum 

In the semiclassical and local approximations, the quantum system corresponding 
to the 2-D anharmonic oscillator introduced in the previous Section is described by an 
effective Lagrangian (effective action density) 



(5.1) 



ie(q,q) == -Z,j{q)q,qj - T4(q), 



where q = ((71,(72) and q = ((71,(72). In the rescaled variables used in the classical case, 
we have 



(5.2) 
(5.3) 

(5.4) 

(5.5) 



Ve = \i<il + 92) + qfqi + 1 ( v/A+ + v/A- - 2 ; 



zu^i+.ri 



{i + v? , {i-v? 



A 



5/2 



A 



5/2 



+ Sqk 



Z12 — Z21 — 7 ■ 



91 92 



(i + 7?)(i + e) , {i-v)(i-o 



A 



5/2 



5/2 



^22 = l+7<^f 



(1+0' , (1-?)' 



A 



5/2 



A 



5/2 



A 



+ 89^C } , 



+ 8qiq2C 
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where 

(5.6) K± = l + ql + ql±Ts, 



(5.7) E = ^(7f + g| + 14g2g|, 

(5.8) ^^{ql + Iql)/^, 

(5.9) e - (92 + 7g?)/S, 
and 



(5.10) 



[(9?+9l)/S]^ 



[,/K^KZ{^K^+,/kZ) 



With respect to the classical system, we have an additional parameter 7 = hg/m'^uj^ 
which arises from rescaling h. It can be seen that, when q varies, the effective potential 
and the symmetric kinetic matrix Zij can be singular or complex-valued, unless q is 
constrained inside a certain region. If we limit ourselves to the region where the effective 
potential, in the considered approximation, is convex, then Lg is well defined. 
The rescaled equations of motion corresponding to the Lagrangian (5.1) 

(5.11) Ziiqi + Zi2q2 = -t;^ — 91 91 7^ — 9i92 + t;^ t: — 9292 - t;— 

dq2 J oqi 



\iqi + Zi2q2 = 


1 dZn . . 

o a 9i9i - 
2 dqi 


dZn . . , 1 

a 91 92 + 
dq2 


( 1 5^22 
^2 dqi 


1291 + ■^2292 = 


1 9^22 . . 
-7:^ — 9292 - 

2 992 


dZ22 . . 1 

- -^ — 9192 + 
dqi 


(IdZii 
^2 dq2 



,r ^n\ r, ■■ , ry ■■ -"".i^. . ""ZZ. . , 1 -""11 dZl2\ . . dVg 

(5.12) Zi2qi + Z2292 = -77-^ — 9292 - -^ — 91 92 + 7;-^ ^ — 9i9i - -^— 

2 dq2 dqi \2 dq2 dqi ) dq2 

have been numerically solved as in the classical case. These equations arc nonlinear and 
may lead to a chaotic evolution. However, due to the fact that no chaotic behavior is 
allowed at quantum level, we expect a reduction of chaoticity in the effective system with 
respect to the classical one. This reduction should depend on the value of the parameter 
7, the value 7 = corresponding to the classical system. In Fig. 1 we illustrate, for 
different values of 7, the smallest energy (threshold energy) eth at which chaos shows 
up in the Poincare sections of the effective system [30]. We see that eth increases with 
increasing 7. This behavior can be explained as follows. Let us consider the Taylor 
expansion of the effective Lagrangian around q = and q = 0. Up to quadratic terms, 
we obtain Le(q, q) = 5(91 +92)^ i(l + 7)(9i +92) + • • •• In terms of unrescaled variables 
this corresponds to a shift of the classical frequency uj -^ lo\/\ + 7. The rescaled energy 
£ — E / {w? u)'^ / g) picks up a factor (1 + 7)^. This means that if et/i(0) denotes the 
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Fig. 1. - Threshold energy eth at which chaotic behavior shows up in the Poincare sections of 
the effective system versus 7 = hg/m^uj^. The solid line is the theoretical estimate £4(1(7) — 
Sth{0){l + 7)^. The dashed line represents the maximal energy Sm below which the effective 
potential is everywhere convex. 



threshold energy at 7 = 0, we should have approximately £4/1(7) = et/i(0)(l +7)^. This 
prediction is well confirmed in Fig. 1 . It is easy to see that the increase of the threshold 
for chaos holds under the general condition that T4i is convex, which, in turn, amounts 
to 



(5.13) He(tr(He(y))^j > 0. 

where V is the classical potential. These results parallel those obtained in [31] for the 
iV-component </>'* oscillators where the mean field plays the role of classical system. 

The range of 7 values explored in Fig. 1 includes situations encountered in molecular 
physics. In fact, the vibrational Hamiltonian of diatomic molecules is often assumed as a 
quartic oscillator and using the numerical values of Ref. [32] obtained from spectroscopic 
data we find that 10"'' < 7 < 10"^ 

In the following we compare the solutions of the local effective equations with the 
classical solutions and with the exact quantum evolutions of coherent states centered at 
the initial conditions of the local equations. We have already remarked that the initial 
conditions for the classical and the effective dynamics do not determine completely the 
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initial wave function but provide only a constraint. Therefore the choice of the initial 
wave function is not unique. A natural choice is represented by a harmonic coherent state 
which is parametrized by the expectation value of position and momentum. In fact, by 
performing simulations with initial wave functions which satisfy the proper constraints 
but arc of arbitrary shape we find that the agreement between the effective and quantum 
dynamics is very poor when the shape of the initial wave function differs substantially 
from that of a coherent state. 

In rescalcd units, the exact quantum dynamics is defined by the Schrodinger equation 

d 

(5.14) i-f—'ip{qi,q2,t) = Htl;{qi,q2,t), 

with 

(5.15) H = ^ipj+Pl) + l{q! + ql) + qlcil 

where pj — —ij-^ and qj = qj, for j = 1,2, are the rescaled momentum and position 
operators. In order to solve (5.14) we represent the rescaled Hamiltonian operator (5.15) 
in the basis of the eigenstates of the associated 2-D harmonic oscillator Hq = \{p\ + 
P2) + 5(91 +92)- The corresponding infinite matrix is truncated and then diagonalized 
with standard techniques [23] . As initial state we choose the coherent state 



(5.16) Ip'i'?1p2'?2> = e-^«iPie^Pi«ie-^«^P^e^P^'J^|Oo), 

where |0o) is the ground state of H{). The parameters p[q'iP2'l'2 ^-re taken equal to the 
initial conditions used in the integration of the classical and effective Lagrangians. 

In Fig. 2 we show the evolution of {qi(t)) in comparison with the corresponding 
classical and effective solutions. The three panels correspond, from top to bottom, to 
increasing values of 7 at constant dimensionless classical energies e. Note that the initial 
coherent state depends on 7. 

Figure 2 shows that by increasing 7 there is a crossover in the behavior of the solution 
of the effective dynamics. At small 7 the effective solution stays close to the classical 
one while for larger 7 it reproduces qualitatively the shape of the quantum evolution. 
We notice that in the large-7 region the quantum and the effective dynamics do not 
show, on the time scale considered, a transfer of energy among the degrees of freedom 
as the classical solution. This seems to indicate that the quantum corrections in the 
effective dynamics have an anti-mixing influence that regularizes the motion. Of course, 
over longer times a transfer of energy takes place also in the quantum and effective 
evolutions. The theoretical implications of these results will be discussed elsewhere. 
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Fig. 2. - Time evolution of the expectation value of the position operator {qi{t)) (solid line) 
compared with the classical (shaded area) and the effective (dots) solutions gi(t). The three 
panels correspond, from top to bottom, to 7 = 0.01, 0.1, and 1, respectively. In all cases we 
have a classical rescaled energy e = 0.1. 
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